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In one-component abelian sandpile models, the toppling probabilities are independent quantities. 
This is not the case in multi-component models. The condition of associativity of the underlying 
abelian algebras impose nonlinear relations among the toppling probabilities. These relations are 
derived for the case of two-component quadratic abelian algebras. We show that abelian sandpile 
models with two conservation laws have only trivial avalanches. 
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I. INTRODUCTION 

Sandpile models are important toy models to under- 
stand self-organized criticality [l| . In an abelian sandpile 
model, the toppling rules can be encoded in an abelian 
algebra. This was pointed out by Dhar0, [||, H| • The 
structure of these algebras [f| (see also is very sim- 
ple, their physical relevance will be shown later in the 
text. They are defined by taking graphs (see Fig. 1) and 
attaching to each vertex a generator <n of the algebra. All 
the generators commute with each other. Two vertices 
are connected by at most two links oriented in opposite 
directions. To each vertex "i" we attribute a polynomial 
relation which expresses a power of (say n) as a poly- 
nomial in the generators attached to the sites reached by 
the outgoing arrows starting at " i" as well as a, . The de- 
gree of the polynomial is at most equal to n. This implies 
for example that for the vertex "0" in Fig. 1 we have 



In the corresponding sandpile model, a§ is interpreted 
as having k grains of sand on the vertex 0. The coeffi- 
cients in the polynomial P(oq, 0,1,0,2,0,3, 04) are nonneg- 
ative and their sum is equal to 1, as we are going to see, 
they are going to be interpreted as probabilities. 

A simple and very relevant example [l[ is the case of 
a two-dimensional square lattice (coordinates the 
abelian algebra being given by the relations: 



a i,j + l Oi— lj a i.j-l, [f 



0. (2) 



We didn't specify the boundary conditions. 

The sandpile model is defined by a stochastic process 
which gives the stationary state and the rules how the 
sand grains act. 

In continuous time and a lattice with N vertices (sites) 
the time evolution of the system is given by a Hamilto- 
nian H : 



N 



P(a , 0,1,0,2,0,3,0,4). 



(1) 



h = y]wj(i - a,i) , 



(3) 




FIG. 1: A typical graph on which one can define abelian 
algebras useful for sandpile models. There are 4 oriented links 
leaving the vertex 0. The sites 2 and 3 are connected by two 
edges of opposite orientation. 
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where the nonnegativc coefficients Wi are transition rates 
and we chose the unit of time by taking: 



JV 

E 



W; 



(4) 



The Hamiltonian §3§ acts in an il/-dimcnsional vector 
space given by all the independent M monomials in the 
generators a^. There is a correspondence between each 
monomial and a configuration in the sandpile model in 
which to each generator appearing in the monomial corre- 
sponds sand grains on the respective sites. Their number 
equals the power of which the generator appears in the 
monomial 

The action of H on the configuration space can be 
understood in the following way: take a configuration, in 
a unit of time, with a probability wi a grain of sand is 



2 



neighbor interactions only are 



FIG. 2: A one-dimensional acyclic directed lattice with N 
sites. All edges have the same orientation. 



added at the site " i" . As a result this configuration goes 
to others, with probabilities given by the toppling rules 
©■ 

The unnormalizcd probabilities P m (t) to find the sys- 
tem in a configuration " to" at the time t can be obtained 
from the master equation: 



dP(t) 
dt 



-HP(t) 



(5) 



where 



M 



P(t) = J2 p m(t)W(m), 



(0) 



W(m) is the monomial corresponding to the configura- 
tion " to" . 

The stationary PDF will be denoted by |0) (H\0) = 0). 

Except for positivity there are no supplementary con- 
straints in the relations ([T]). This implies that one can 
choose at will the toppling probabilities, however as we 
are going to see this is not the case if one considers two- 
component sandpiles models. 



II. SANDPILE MODELS WITH TWO KINDS OF 
SAND 



In two-component abclian sandpile models, one as- 
sumes that one has two types of sand (say "a" and "6"). 
One adds with a probability Ui a grain of sand "a" 
(respectively " b" ) to the site " i" . Grains topple if there 
are more than one, of any kind, on a given site. In the 
toppling process, the two types of sand mix and eventu- 
ally transform into each other. We are going to consider 
the simplest model of this kind. 

Consider a one-dimensional directed lattice (Fig. 2) 
with N sites To each site "i" we attach two gen- 
erators a, and 6^, all of them being mutually commuting 



\Oi, a v 



[a h b v ] = KM = Vi,i' = l,...,N. (7) 



For simplicity we consider quadratic algebras only. 
The most general quadratic relations involving nearest- 



a\ = ai a t a l+1 + a 2 a t b t+1 + fli bia i+ i + (3 2 b t b t+1 
b\ = 7i a t a l+1 + 72 a;6;+i + Si bia i+1 + S 2 hb i+ i 

+ vi a i+i + m b 'i+i + vsai+ibi+i, 

ctibi = hi Oifli+i + fj,2 a,ib l+ i + vi bia l+ i + v 2 b,b l+ i 
+ Cia?+i + (2b 2 l+1 + C, 3 a i+ ib i+ i, 

(8) 

for i = 1, . . . , N, in which we take: 



a-N+i — bjy + i — 1 



(9) 



This implies that the two types of sand may leave the 
system on the site N. The constants in © arc positive 
probabilities, so that 

« + + £ = J + S + rj = Ai + ^ + C = 1, (10) 

where we used the notation a = Y^i=i a i> P = Ei=i Pi: 

6 = Ei=i&, ctc - 

The relations ([T])-© don't define yet an algebra since 
one has still to impose associativity (diamond condi- 
tions): 

(of)6i = a,i(aibi) , a,(6?) = (a i 6 i )6 l . (11) 
Introducing (jSJ in (jlip we find the following 12 relations: 



Ci = Pi{m - 5i) + Vi{i>i -oti), i = l, 2, 

6 = Pl(p2 - h) + faipi - 81) 

+vi(v2 - a 2 ) + v 2 (vi - ai) , 

m = li{vi - oti) + fXi{fJbi - Si) , i = 1, 2, 
m = 71(^2 - aa) + 72(^1 - oli) 

+Mi(M2 - S 2 ) + a*2(mi - Si) , 

C» = fali - , i = 1,2, 

Cs = /?i72 + /?27i - Mi^2 - ^2^1 , 



(12) 



and 



13(1+ n-S) = (l-v)(l + v-a), (13) 
7(1 + 1/ -a) = (1-^(1 + ^-5), (14) 
/3 7 = (15) 



Out of the 3 relations p3|) - p5| only two are independent 
since multiplying (flU)) with (fT4")) one obtains (fT5|) . 

As one can see from (fT2" |) - (fT5"|) in the two-component 
case the various toppling probabilities are constrained by 
nonlinear relations. 

The vector space in which the Hamiltonian 



H 



1 N 

-y J~] (1 - Ui Oj - Vi bj) , Ui + Vi = l, (16) 



4=1 



acts is made out of monomials in which cither a j , 6j or 1 
appears for a given site " i" . The physical interpretation 



3 



of these monomials is obvious. If in the monomial at (bt) 
appears, on the site "i" one has a grain of sand of type 
"a" ("6"). If neither or hi appear in the monomial, 
there is a vacancy on the site " i" . 

It is easy to show the the stationary state PDF is of 
product form: 



N (1- 



a)(l-5)-fj.i>+(l+fj.-S)ai + (l+v-a)bi 
(2-a)(2-S)-(l- M )(l-i>) 



(17) 



io) = nr=i 



An obvious question is if there is a solution of (|I2 ]) -(|Io |) 
in which each type of sand is conserved separately. This 
would imply perhaps a new universality class of sandpile 
models. If one takes only a\, £i, 82, T72, fJ*2> v i an d C3 as 
non vanishing probabilities, in the bulk, the number of 
sand grains of type " a" and " b" are conserved separately 
There are only two equivalent solutions of the diamond 
conditions in this case. One solution is: 

a\ = aa,» H i + (1 - a) a? +1 , (18) 
bj = bib i+1 , (19) 
CLibi = b t a l+1 . (20) 

In the second solution, one exchanges the generators a, 
with 6j. 

The stationary PDF of the stochastic process (JTBJ), 

(nsj-jnii is 

N 

\0)=l[b i . (21) 
»=i 

There are only trivial avalanches in this case. If a grain of 
sand of type " a" hits the system, it leaves directly at the 
boundary because of (f^Uj) . leaving the system unchanged. 
If a sand grain of type " 6" hits the system, it also leaves 
it directly at the boundary because of ([19]). This is a 
surprising result. 

We consider now a more general graph where the ver- 
tex " z" is linked by outgoing arrows to a number of sites 
which wc label by index " x" . If one asks for two conser- 
vation laws and imposes the associativity condition (|lip , 
the generalization of the solution psp - pp]) reads: 

a* = 22 a * a i a x + ^2^x, v a, x a y , (22) 

x x : y 

bj = ^Svbib*, (23) 

X 

aibi = 2J v x ha x , (24) 

X 

where the toppling probabilities satisfy the relations 
J2 X $x = XL v * = 1 ' a * - v * ' and 



; ■= 5 W 1 ^ - OL v ) + v y (v x - a x )}. (25) 



Notice that, similar to the one-dimensional case, only 
the grains of sand of type "a" topple. This has as conse- 
quence that for an arbitrary graph the stationary state of 



the system is an absorbing state, like (|2~lj) , in which only 
particles "b" are present. The proof is straightforward. 
Only trivial avalanches can be obtained. 

We have shown that using other representations of the 
algebra ([T8|) - (j20|) which may define other sandpile models 
(see 0), one obtains the same result: two conservation 
laws are compatible with trivial avalanches only. The 
system stays unchanged in the aftermath. 

We have checked, only for quadratic algebras, that for 
a vertex connected by arrows to several vertices (there 
was only one in example ©), the associativity condi- 
tions (|lip are incompatible with the existence of non- 
trivial avalanches in the sandpiles with two conservation 
laws. This is a surprising result. We have no proof of a 
similar statement for more general algebras (cubic, quar- 
tic,. . . ). Simple examples like the generalization of 
give the same result. 

If we don't insist on two conservation laws, there are 
various solutions of the constraints (|T2|) — (fT5j) and there- 
fore two-component sandpile models. They all share the 
property that during the toppling process, grains of sand 
of type "a" and "6" mutate (we have not identified a 
physical process which leads to such a phenomenon). 
We present a simple example in which we consider the 
abelian algebra: 



a] = (1 - v) a t a l+ i + 



1 + /i 



bidi+i 



+ TT^+ 1 + ^T7r ai + l6 '+ 1 ' (26) 

b\ = (1 - pl) kb i+ i + 0(1 + (J, - 2(f)) a i+1 b i+ i 

+ 2 a? +1 + (l-0)(/i-0)^ +1 , (27) 

a^i = 4>ai(ii + i + (1 — </>) ai&i+i . (28) 



and the Hamiltonian (|16p with = = 1/2. This 
implies that we take the stationary state 



N 



i°>=n 



fiv + (1 + n)<Xi + vbi 
(1 + + 



(29) 



and add with equal probability a grain of sand of type " a" 
and "6" on the first site. Notice that |0) is independent 
on 0. 

We are interested to know what is the probability 
Pa (T) to have an avalanche ending with a grain " a" and 
having a duration T. The duration of the avalanche, in 
this one dimensional case, is the number of sites where 
the topplings occur. P&(T) has a similar meaning. In Q 
it was shown that in the case of the one-component model 
which is obtained by considering only (|27p in which one 
takes </> = 0, at large values of T, one has 



P b (T) - T~ 3/2 . 



(30) 



This means that for the one-component case the 
avalanches are in the random walker universality class 
Q. In Fig. 3 we present the results of Monte-Carlo sim- 
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FIG. 3: (Color online) Probability distribution functions 
P a {T) and Pt{T) to have avalanches with duration T end- 
ing with grains of sand of type "a", respectively, "6" obtained 
in the model given by (|26[) - (|28[1 , The values [i = v = 1 and 
4> = 1/2 were used for the parameters of the model. 4 x 10 9 
avalanches were observed in the Monte Carlo simulations. In 
the figure P a (T) and P b (T) are multiplied by T 3/2 in order to 
show their large T behavior. 



ulations for the model given by (|2"6")) - (f2"5)) in which we 
have taken fi = u = 1, 0=1/2. We observe that for 
large values of T. 



P a (T)~T~ 3/2 , P b (T)~T- 



3/2 



(31) 



toppling probabilities are not arbitrary. They have to 
satisfy nonlinear relations coming from the condition that 
the algebra is associative. These constraints don't exist 
in the one-component models. An unexpected result is 
that, at least for the case of quadratic algebras, the co- 
existence of two conservation laws in the bulk (one for 
each component), and of nontrivial avalanches, is impos- 
sible. 

We have also shown, in an example of a one- 
dimensional directed model, that once we allow the two 
components to mutate in each other during the toppling 
process, the one-component and the two-component 
models belong to the same universality class. 

We believe that our conclusions apply to any multi- 
component models. 

We would like to mention a possible extension of the 
quadratic algebra (8). If we omit the last of the three 
equations (8), two grains of sand, one of type "a" and 
the other one of type "b", don't topple. As a result the 
stationary state is a linear combination of 4^ instead of 
3 N states. If we are not interested in mutations (a <-> b), 
one has a direct product of two algebras. One contains 
only dj generators (like (18)) and the other one contains 
only hi generators. One gets avalanches in which grains 
of type "a" and "b" don't mix. If one consider the pos- 
sibility of mutations the situation is different. Take for 
example the algebra defined by (26) and (27) (we have 
omitted (28)!), the total number of "grains" is not con- 
served since a^+i^i+i has to be looked upon as a new kind 
grain. As a result, the PDFs of the duration of avalanches 
ending in a grain of sand of type " a" , " b" or " ab" have 
an exponential falloff. We have also checked that if we 
change the algebra by allowing mutations, while conserv- 
ing the number of sand grains, the avalanches belong to 
the random walker universality class 0] • 



mutations between sand grains of type " a" and " b" don't 
change the universality class. This stays valid for the 
whole parameter space except for the boundaries = 
and <f> = /.t. 

If in ip5 ]) -([28"] ) one takes = 0, P a (T) decreases ex- 
ponentially and Pb{T) has an algebraic falloff (|50|) . This 
can be easily shown using the algebra (|26"|) - (f28]) and the 
recurrence relations for the toppling probabilities which 
follow. Conversely, if = fx, P\(T) decays exponentially 
and P a (T) has the algebraic falloff. 



III. CONCLUSIONS 

The main message of this paper is that if one is in- 
terested in two-component abelian sandpile models, the 
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